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Resumen 

En este trabajo, se obtiene una expresion numerica para 
calcular el espectro de emision de un sistema microcavidad- 
punto cuantico usando una teon'a de campo medio en el 
formalismo de la matriz densidad. El sistema modelado es 
un micropilar semiconductor que contiene un linico punto 
cuantico en el interior de la microcavidad, este sistema pre- 
senta perdida de fotones a traves de los espejos de la cavidad 
y bombeo de excitones al punto cuantico. Obtenemos una 
ecuacion maestra de campo medio no lineal que nos permite 
calcular el espectro de emision. 
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Abstract 

In this work, we deduce a numerical expression to calculate 
the emission spectrum for a quantum dot-microcavity sys¬ 
tem, using a mean field theory in the density matrix forma¬ 
lism. The open system modeled, is a typical semiconductor 
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micropillar that contains a single quantum dot inside a mi¬ 
crocavity, this system presents leakage of photons through 
cavity mirrors and pumping of excitons to the quantum dot. 
We obtain a mean field nonlinear master equation that allow 
us calculate the emission spectrum. 


Keywords: Microcavity, quantum dot, mean field theory 


Introduccion 

Desde que Purcell descubrio el efecto que una cavidad tiene en el 
tiempo de relajacion de un emisor optico [1], lo emisores dentro de 
cavidades han sido de gran interes investigative debido a sus po- 
tenciales aplicaciones tecnologicas, especialmente desde su imple- 
mentacion en sistemas de estado solido como lo son, en particular, 
puntos cuanticos al interior de micropilares. 

Es bien conocido que en estos sistemas se presentan dos regimenes: 
acople debil y acople fuerte (WC y SC, respectivamente, por sus 
siglas en ingles). En el primero, la interaccion entre la radiacion 
(modes opticos de la cavidad) y la materia (puntos cuanticos) es 
tan debil que puede ser tratada perturbativamente; en el ultimo se 
presenta una alta probabilidad de emision y reabsorcion de fotones 
por parte del emisor, generando estados altamente acoplados de luz 
y materia conocidos como estados vestidos. |2] 

Un modelo usualmente utilizado para la descripcion de la fisica 
fundamental de este sistema es el de Jaynes Cummings {h = 1)[3] 

H = UfO^a + UaO'^o'+ g{a^a + a^a) (1) 

a y son los operadores de aniquilacion y creacion de fotones de 
la cavidad, respectivamente, con energia ojj, estos siguen la regia 
de conmutacion [a,a^ = 1; a y son los operadores escalera para 
exciton, ademas satisfacen la regia de anticonmutacion {a, = 1 

debido a la estadistica fermionica que le imponemos al punto cuanti- 
co de energia oJa] g indica la fuerza de acople lineal entre radiacion 
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y materia. 

El sistema que vamos a tratar incluye, ademas, una perdida de fo- 
tones a traves de los espejos de la cavidad, una perdida de excitones 
debida a emision espontanea y un bombeo incoherente de excitones 
al punto cuantico. Nos interesa analizar la validez de la aproxima- 
cion de campo medio en la cual el valor esperado del producto de 
operadores se expresa como el producto de los valores esperados de 
cada operador[lj. Luego de analizar la validez de la aproximacion 
se propone un calculo del espectro de emision. 

Fijamos g = ImeV y lo tomamos como nuestra escala de energia[S] 
a partir de la cual se determinan las demas cantidades. 


Marco Teorico 

Estudios previos ban mostrado la forma de tratar este sistema 
afiadiendo terminos disipativos de Lindblad a la ecuacion maestra 
de Lionville-von Neumann para la evolucion temporal de la matriz 
de densidad[5][6] 


dtp = i[p, H] + ^(2apa^ — a^ap — pa^a) + ■^{2apa^ — a^ap — pa^a) 

p 

+—{2a^ pa — aa^ p — paa^). ( 2 ) 


7 , 


El sistema pierde fotones con una tasa k a traves de los espejos de la 
microcavidad, 7 es la tasa de decaimiento del exciton por emision 
espontanea y P es la tasa de bombeo continuo e incoherente de 
exciton al punto cuantico. Calculando los elementos de matriz en 
la base {\Xn)] \Gn)) se obtiene pin,jm = lo cual lleva al 

siguiente conjunto infinite de ecuaciones diferenciales[5]: 


dtPGn,Gn —*5'\/^(PGn,Xn-l “ PXn-l,Gn) + 'lPXn,Xn 

l^\jlpGn,Gn P l)pGn+l,Gn+l] PpGn,Gn 


( 3 ) 
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dtPXn,Xn —igVn + l{pXn,Gn+l — PGn+l,Xn) — lPXn,Xn 

— K[npxn,Xn — {n + l)pxn+l,Xn+-\\ + P PGn,Gn (4) 


dtpGn,Xn-l —i[gVn{pGn,Gn “ PXn-l,Xn-l) “ ^pGn,Xn-l] 

- [(7 + K{2n - 1) + P)/2]pGn,Xn-l 

- K^/n{n + l)pGn+l,Xn (5) 


dtPXn-l,Gn — — i[g\/n{pGn,Gn “ PXn-l,Xn-l) “ ^PXn-l,G'n\ 

— [(7 + K(2n — 1) + P) /2]pxn-l,Gn 

- K7/n(nTI)PXn,Gn+l (6) 

donde se ha definido el detuning A = Uf—uja- Una vez calculados los 
elementos de matriz, el niimero medio de fotones se puede obtener 
de: 


N = ^ rn{pGm,Gm + PXm,Xm) (7) 

m 

Este tratamiento es lo que consideramos como ’’modelo exacto”, 
ya que su planteamiento es anahtico a pesar de que sea necesario 
truncar el conjunto de ecuaciones diferenciales; respecto a el com- 
pararemos nuestros resultados. 

Teniendo en cuenta que el valor esperado de un operador se obtie- 
ne como (O) = Tr{pO) y la derivada del mismo en el cuadro de 
Schrodinger es (O) = Tr{Odtp), se calcula la dinamica de los valo- 
res esperados haciendo uso de ([^. A1 hacer esto, se obtienen valores 
esperados para los conjuntos de operadores (crVa'^) y (a^aa^a), los 
cuales en la aproximacion de campo medio se expresan como: 


(a^aa^) ~ {a^a){a'^) 


( 8 ) 
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(a^acrV) ~ {a^){a){a^a) 

( 9 ) 

obteniendo las siguientes ecuaciones: 


(at) = iujf{a)) + ig{a^) - ^(at) 

( 10 ) 

(at) = iu:a((y^) + ig{a)) - 2i5f(at)(ata) - ^(at 

) (11) 

(crta) = ig{{a)){a) - (at)(a)) - r(ata) -h P 

( 12 ) 

y los hermiticos conjugados de (10) y Se ha definido V = 

7 -|- F una tasa de disipacion efectiva para el exciton. La razon para 
realizar la aproximacion asi y no, por ejemplo, como (ataat) cs 


{a es porque de esta forma no se pierde informacion del 

bombeo P el cual es importante para mantener las excitaciones en 
la cavidad, permitiendo la llegada a un estado estacionario, veamos. 
Si se tomara (aV) ~ (a'*')(a), su derivada seria (crta) ~ (crt)(cr) + 
(a'*')(a), con lo que se obtendria: 

(atcj) ig{{a)){a) - (a^)(a)) - r(cT^)(cT) (13) 

Comparando con ( [I^ vemos que al separar (crV) se pierde parte 
del efecto que tiene el bombeo sobre el sistema. 

Con la dinamica de los valores esperados resuelta, se define la fun- 
cion de correlacion en la aproximacion de campo medio, en analogia 
a lo realizado en [2], como: 

G'“^((,r) = (a'(«))(a(« + r)> (14) 

Con la funcion de correlacion se puede calcular el espectro, en par¬ 
ticular, se puede calcular para estados estacionarios que son en los 
que nos vamos a centrar: 

Siuj) = — lim TZ 

M t^oo 

siendo M una constante de normalizacion. 


G"''(t,T)e“^dT 


( 16 ) 
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0.20 

/{•(meV) 



Figura 1. Numero medio de fotones en aproximacion de campo medio N en 
escala de grises, como una funcion de k y P en condiciones de resonancia 

A = 0 con 7 = 0.1 y g = 1. 


Validez de la aproximacion 

A1 realizar la aproximacion se pierde informacion de las reglas de 
conmutacion debido a que al proponer (a'^a) ~ (a'*')(a), los valo- 
res esperados del lado dereclio conmutan, equivalente a decir que 
[a, a)] = 0. Es de esperarse, entonces, que la aproximacion dependa 
de las poblaciones medias de fotones en la cavidad, lo que la liace 
sensible a los parametros k j P. 

Definimos el numero medio de fotones en la aproximacion de campo 
medio como: 


N = {a^){a) 


( 16 ) 
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Figura 2. Error relativo 5, como una funcion de k y P en condiciones de 
resonancia A = 0 con 7 = 0.1 y g = \, se observa un bajo error para pequenos 

valores de disipacion. 


La figura muestra como varfa este ultimo en funcion de kj P. Se 
obtiene un alto mimero de fotones para K?«0.1yP~20, algo que 
es de esperarse ya que este disminuye con su tasa de disipacion y 
aumenta con el bombeo de excitaciones en la cavidad. Sin embargo, 
el bombeo del sistema se realiza en el exciton, de forma que valores 
de P muy altos producen una saturacion en la poblacion de excito- 
nes, lo que impide el crecimiento de la poblacion de fotones como lo 
muestra la figura, comportamiento que ya ha sido estudiado antes 

01711 . 

Con el objetivo de realizar un analisis cuantitativo de la validez de la 
aproximacion de campo medio y de su sensibilidad a los parametros 
disipativos, se define el error relativo 6 = |(iV — iV)/iV|. La figura 
[^evidencia el comportamiento de esta cantidad en funcion de k y 
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GURA 3. 


Espectro de emision del sistema en resonancia para los pardmetros 
K = 0.1, P = 10, 7 = 0.1 y g = 1. 



a;(meV) 


Figura 4. Espectro de emision del sistema con detuning A = —5, n = 0.1, 

P = 10, 7 = 0.1 2/ 5 = 1. 
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P. Se encuentra una region bien definida en la que el error relativo 
es menor del 10 %. Vemos, ademas, que para grandes valores de 
K, 6 crece, esto es debido al bajo niimero medio de fotones que se 
presenta. Como estamos olvidando las reglas de conmutacion, es de 
esperarse que la aproximacion tenga mayor validez a medida que 
aumenta la poblacion de fotones en la cavidad. 


Espectro de Emision 


Ya conociendo las regiones en que la aproximacion es valida, nos 
remontamos al calculo del espectro de emision, que consiste en la 


transformada de Fourier de la funcion de correlacion (14). El es¬ 


pectro obtenido, figura presenta su maximo pico en frecuencia 
u ~ lOOOmeV, algo que se esperaba debido a los valores tornados 
para las energias propias de la cavidad ojf^a- La figura muestra el 
efecto que tiene el detuning. Los espectros fueron calculados para 
estado estacionario. 


Conclusiones 

Se analizo cualitativa y cuantitativamente la aproximacion de 
campo medio aqui descrita, para el calculo del mimero medio de 
fotones y el espectro de emision de un sistema microcavidad-punto 
cuantico. Se encontro que la aproximacion es valida para valores de 
K entre OmeV y 0.4meV, y para P entre ImeV y 14meV, con un 
error relativo inferior al 10 %. Se calculo el espectro y este reproduce 
el resultado conocido de reduccion del ancho de linea para el caso 
de alto bombeo de excitones en la cavidad. 
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